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ABSTRACT 



^ , A search for possible radiative decays of solar neutrinos with emission of photons in the 

O ■ visible range was performed during the total solar eclipse of August 11, 1999. Due to very 

■ bad weather conditions our two telescopes were unable to collect useful data; fortunately we 

^ . obtained several video camera images from a local TV station. An analysis of the digitised 

Qh! images is presented and limits on the lifetime for radiative decay are discussed. 

1. Introduction 

^1 It is generally agreed that most probably neutrinos have non-zero masses. This belief is 

' based primarely on the evidence/indication for neutrino oscillations from data on solar and 
atmospheric neutrinos. 

Neutrino oscillations are possible if the flavour eigenstates are not pure mass eigenstates, 
e.g.: 

\ve >= Wi > COS 9 + |z/2 > sin 6* (1) 

where rriiy^ > ^i^i ^i^d 6 is the mixing angle. 

Since few years there is evidence that the number of solar neutrinos arriving on Earth 
is considerably smaller than what is expected on the basis of the "Standard Solar Model" 
and of the "Standard Model" of particle physics, where neutrinos are massless (see e.g. [jl|). 
One possible explanation of these experimental results involves neutrino oscillations, either in 
vacuum with Amf^^ = ml^ — ~ 10~^° eV^ (as originally discussed in [0, |^) or resonant 
matter oscillations Amf^^^ ~ 10"^ eV^ 0, ||. 

Recent results from Super-Kamiokande 0, MACRO [0] and Soudan2 |^ experiments on 
atmospheric neutrinos support the hypothesis of neutrino oscillations (in particular — > v/) 
with large mixing (sin^ 29 > 0.8) and Am^^„ ^ 3 x 10~^ eV^. 
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Figure 1 : Sketch of the principle of the experiment to detect radiative neutrino decays during a solar 
eclipse (the emission angles of the photon and of the neutrino are enlarged). 



Another indication in favor of neutrino oscillations with a third energy scale Aml^^^, ~ 1 
eV^ is reported in 0. 

The above observations appear to be the first indications for new physics beyond the 
"Standard Model", and any model that generates neutrino masses must contain a natural 
mechanism that explains their values and the relation to the masses of their corresponding 
charged leptons. 

Different scenarios have been proposed to explain all the observations, including the results 
with neutrinos from reactors and accelerators |T^. 

If neutrinos do have masses, then the heavier neutrinos could decay into the lighter ones. 
For neutrinos with masses of few eV the only decay modes kinematically allowed are radiative 
decays of the type Ui ^ Uj + •y (where lepton flavour would be violated). 

Upper bounds on the lifetimes of such decays are based on astrophysical non-observation 
of the final state 7 rays. Limits were obtained from measurements of X and 7 ray fluxes from 



the Sun [131 and SN 1987A M 



In the case of neutrinos with nearly degenerated masses, of the order of the eV, the emitted 



photon can be in the visible or ultraviolet bands [0, |T8[. A first tentative to detect such 
photons, using the Sun as a source, was made during the total solar eclipse of October 24, 1995 



Direct visible photons from the Sun come at a rate of some 10^^ cm~^ s^^; this makes a 
direct search for photons from solar neutrino decays impossible. To perform a measurement 
one must take advantage of a total solar eclipse, which cuts by at least 8 orders of magnitude 
the direct photon flux. By looking with a telescope at the dark disk of the Moon one can search 
for photons emitted by neutrinos decaying during their 380000 km flight path from the Moon 
to the Earth, Fig. |I|. 

In this paper we describe the methodology employed for such a search and give limits 
obtained from a preliminary measurement performed during the total solar eclipse of 11 August, 
1999. 

2. Kinematics of radiative decays 

We assume the existence of a possible neutrino radiative decay, 1^2 ^^1+7, where rriu^ > rriy^ 
and z/i, z/2 are neutrino mass eigenstates. 

The energy of the emitted photon in the earth reference laboratory system is 
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Figure 2: (a) Expected maximum angle of emission of visible photons from radiative solar neutrino 
decays as function of the vi mass, (b) Percentage of visible photons for different z/i mass values. 
The open triangles correspond to the Super Kamiokande solution with Am^ = 6 x 10~^ eV^ and 
sm^26 = 4 X lO^'^; the black points correspond to the solution with Am? = 2 x 10~^ eV^ and 
sin^ 29 = 0.71. 



Elab = Ecmlu (1 + COS 6*] 



where Ei, and 'j^ = = (1 — /?^) ^ are in the lab. frame, and 6* and E^m are the photon 
emission angle and the energy of the emitted photon in the decaying neutrino rest frame. 
For radiative neutrino decays the general expression for the angular distribution of the 



emitted photons in the rest frame of the parent neutrino is [|T£ 

dN 1 



1-acosr) (3) 
a cos u* 2 

where the a parameter is equal to -1, +1, for left-handed and right-handed Dirac neutrinos, 
respectively, and for Majorana neutrinos. 

In order to estimate the expected number of photons in the visible range and the maximum 
angle of emission of visible photons from radiative solar neutrino decays we performed Monte 
Carlo simulations for neutrino masses in the range 0.1 — 4 eV and two different scenarios of 
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solar neutrino oscillations (Am^ = 2 x 10""^ eV^ and sin^ 2^^ = 0.71; Am^ = 6 x 10"® eV^ and 
sin^ 26 = 3.98 x 10"^). We used the following procedure: 

i) randomly generate the quadrimomentum of 1^2 with reference to the solar neutrino energy 
spectra pO| ; 

ii) randomly generate cos 6'* by using Eq. ^ with a = —1, in agreement with the Standard 
Model; the angular distributions in the lab. (earth) frame were found to be essentially 
insensitive to the choice of the a parameter pl| ; 

iii) associate to each photon a quadrimomentum from the energy-momentum conservation; 

iv) transform the computed values to the laboratory (earth) reference system. 

Fig. H is the result of 5 x 10® simulated radiative decays, for each (mj,^, Am^) set of values. 

3. The experiment 

We implemented an experiment designed to make observations during the total solar eclipse 
of August 11, 1999 (NOTTE, Neutrino Oscillations with Telescope during the Total Eclipse) . 
The aim was to exploit the possible visible photons emitted in 1/2 ^ z^i + 7 decays, during the 
solar neutrino flight from the Moon to the Earth and the shielding of the direct solar light by 
the Moon disk. 

We intended to use two optical telescopes: a 25 cm Newtonian installed in the Parang Massif 
in Romania, close to the point of maximum eclipse; and a smaller one (12.5 cm Cassagrain) 
mounted on an automatic pointing device in a MIG-29 supersonic fighter of the Romanian Air 
Force. The airborne telescope would have had the advantages, compared to the ground one, of 
longer observation times and of reduced hght absorption and diffusion in the atmosphere. 



thus enlarging the detectable photon energy range. 

The telescopes were equipped with CCD's and co-axial digital TV cameras for the control 
of the alignment. The expected overall integrated acceptance (1.3 x 10^ cm^s) could allow 
an improvement of factor at least 20 relative to previous trial [ P^ . 

The extremely bad weather conditions made the observations impossible. We only had the 
possibility to use a SVHS video record offered to us by the television of Ramnicu Valcea; we 
obtained 2750 frames that allowed us to test our procedure for the extraction of a possible 
signal in the digitized frames. 

The SVHS video was obtained with a Panasonic M9000 camera with a focal length of 
200 mm and a diameter of 49 mm. The video was digitized using a 450 MHz Pentium 
n with 256 MByte RAM and a video card Marvel G200 with complex video imput and a 
CoDec (Compression/Decompression) hardware MJPEG which allowed a maximum resolution 
of 704 X 576, with compression of the frames. The used video player was a Panasonic NVHSIOOO 
(standard SVHS) with all the connections in Y/C. The maximum resolution we were able to 
obtain without compressing the frames was 352 x 288 pixels. 

We selected images of 32 x 32 pixels centered in the center of the Moon; the final image 
was the sum of the reduced 2747 available frames (in the following we shall refer to it as the 
measured image). 

4. The analysis procedure 

We used the measured image to test our data reduction software. The measured image is 
dominated by the image of the Moon in the solar light reflected by the Earth. In order to test 
this hypothesis we calculated the linear correlation coefficient between our measured image and 
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Figure 3: Procedure for trying to detect a radiative decay signal in the white channel. We give the 
average luminosity (arbitrary units) as a function of 6 in arcsec for: A, the normalized measured image; 
B, the normalized image of the Moon from the Pises observatory; C, the image after subtraction (= 
measured-Pises); D, the 4*'' order wavelet decomposition image. 



a CCD image of the full Moon (obtained by a group of the Pises observatory [|22|), from which 
we selected the same central region; we used the formula 

VEij {xij - x) Eij ivij - y) 

where the matrix elements Xjj and are the luminosities in each pixel of the two images. 

We obtained a value C{x,y) ~ 0.6, which implies that the two images are positively 
correlated and that we can thus substract one from the other (in the following we shall refer to 
this as the image after subtraction). 

We tested this result by repeating the correlation analysis using the Moon image and 
10^ images randomly generated in accordance with the light distribution in our measured TV 
image; the 10^ images were generated via a Monte Carlo simulation using the HBOOK routine 



HRNDM . We found that the linear correlation coefficient distribution is narrowly centered 
around zero, indicating that the correlation between our measured image and the Pises image 
of the Moon is not accidental. 

We performed an intensity and contrast normalization, by imposing that the pixels average 
intensity and the lowest and the maximum pixel intensity values were the same for the two 
images. This procedure was independently applyed in all (red, green, blue and white) channels. 

On the image after subtraction we made a fourth order wavelet decomposition. Fig. ^ 
illustrates our procedure in the white light; this procedure was performed in all color channels. 
In this way we atenuated any signal component (zodiacal light, corona light, etc. diffused by 
the atmosphere) on a scale greater than 2x2 pixels. 

The expected signal from neutrino radiative decays is concentrated in the central pixel of 
about 20 arcsec. Neutrinos produced through thermonuclear reactions in the inner region of the 
Sun come from a relative small spatial volume corresponding in angular size to about 0.6% of 
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the Sun diameter (~ 12 arcsec); the maximum emission angles for visible photons, with respect 
to the initial neutrino flight direction, would always be limited to few arcsec (see Fig. H). 



5. Estimate of the lifetime sensitivity 
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Figure 4: Average residual light fluxes, after moon subtraction and wavelet decomposition, in the 
red, blue, green channels and summed (white) expressed in photons cm~^ s~^ as a function of the 
angular distance from the center of the Moon disk expressed in arcsec. 

In order to estimate the lifetime of z/2 (see Eq. ^ we need to determine the flux of visible 
photons from z/2 — > i^i + 7 decay. We performed a wavelet decomposition on the image after 
subtraction in all color channels. 

The wavelet analysis is generally used in order to restore the image and to eliminate spurious 
signals. We were interested in what is usually considered as "noise"; thus we did not use the 
terms of wavelet decomposition but we used the residues; in particular we used the fourth 



order residuals of a wavelet decomposition in the Haar basis [^]; the fourth order is the highest 
usable order because the subtracted image has dimensions 32 x 32 pixels Q 

Fig. ^ shows the residuals light signal after wavelet decomposition, in the red, green, 
blue and white channels as a function of the distance from the center of the Moon disk. At the 
beginning of our analysis the 4*'* order residual flux of the wavelet decomposition was expressed 
in Acquisition Digital Units (ADU). We made a conversion from ADU to photon flux units 

"'^We recall that if the original image has 32 x 32 pixels, the 1*'* order image has 2x2 pixels, the 4*'* order 
image has 16 x 16 pixels. 
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using the intensity profile of tlie solar corona. We used as calibration curve the intensity profile 
of the solar corona taken by a team from the Pises Observatory during the same total solar 



eclipse. It presents two landmarks: the intensity of a lunar sea (12100 ADU) and of the Tycho 
crater (38800 ADU). We constructed the intensity profile of the solar corona for our measured 
image of the eclipse of August 11, 1999 in the red, green, blue and white channels as a function 
of the angular distance from the center of the Sun. 

Note that the intensity profiles differ one from the other in the ordinate scale because the 
ADU is a unit which depends on the instrumental characteristics. 

In order to estimate the conversion factor between the two ADU scales we interpolated the 
solar corona intensity profiles, both in our measured image and in the Pises image, with a 
cubic polynomial, and then we integrated them from 1.1 Rq to 1.5 Rq. We obtained a ratio of 
conversion, f=(Pises ADU)/(our ADU)= 8.25. 

The brightness of the full Moon, observed outside the atmosphere, when the Moon is at 
its apogee and the Earth is at its mean distance from the Sun, is B = (0.34 ± 0.01) lux. It 
corresponds to A^^ = (140.0 ±4.4) x 10^ photons cm~^ s~^ at an average wavelength of 5500 A. 
We assigned to this luminosity the intermediate value between the two reported on the Pises 
intensity profile; we estimate an error of about 25% on the conversion factor. 

Fig. H shows that after the moon image removal and the wavelet decomposition of the 
remaining signal, a peak survived in the central pixel only. A conventional optical phenomenon 
that could lead to such a signal is the so called "Poisson spot" (a consequence of the Fresnel 
diffraction) . 

The diffraction occurs if the distance from the light source (the Sun) to an opaque disk (the 
Moon) or from a disk to a screen is much larger than the disk diameter and if the disk is much 
larger than the wavelength; it consists in the appearance of a luminous spot in the centre of 
the shadow (Poisson spot). 

For a point source the intensity of the Poisson spot should be the same as the intensity that 
would be induced by the source at the distance of the screen; for plane waves the intensity of 
the spot is a quarter of that value. In our case it is difficult to estimate the expected intensity 
of the Poisson spot, as the Sun is an extended light source. Considering the light along a line 
perpendicular to an opaque circular disc, we expect to see darkness immediately behind the 
disc; the relative intensity should increase with increasing distance between the disc and the 
point of observation. 

The contribution of the Poisson spot to the peaks in Fig. ^ can be estimated using the 
property that the intensity of the Poisson spot should not depend on the wavelength. Thus the 
ratios between the different Poisson intensities in different color channels should be the same 
as the ratios between the intensity recorded in similar conditions in the same channels of the 
solar spectrum. 

Later we recorded the solar spectrum with a TV-camera similar to the Ramnicu Valcea 
one, and digitized the image in the same way as the measured eclipse image. 

The TV-camera was equipped with a CCD, which had three different physical pixels for 
every pixel of the image: the first pixel is without filter while the others have red and blue 
filters. The green channel is obtained from a subtraction from the white signal of the summed 
red and blue ones; so it is in fact useless for the analysis described below. 

From the solar spectrum frames we obtained a maximum signal of 56 ADU in the red field 
and of 88 ADU in the blue field; both with a background of 14 ADU. Fig. El gives in the central 
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Table 1: Preliminary lifetime lower limits for a possible radiative decay of solar neutrinos, obtained 
from the 1999 total eclipse images. The limits are given as a function of m,yj, Am^ and sin^ 26. In the 
first three columns are given the oscillation parameters used to obtain the limits presented in the last 
two columns. 

pixel a flux = 7.4 x 10^ cm-^ in the red channel and = 1.4 x 10^ cm-^ in the 
blue one. 

We ascribe to the Poisson spot all the signal in the red channel; thus in the blue channel 
we obtain a residual photon flux of about ~ 7% of the signal. 

We assume that the relative contribution to the overall decay signal in the white light may 
be taken as the average of the signals at the two extremities of the visible spectrum, that is 
3.5% of the white maximum, thus = 7.4 x 10^ cm-^ s~^. 

Preliminary results on the limits of the z/2 lifetime, with respect to its radiative decay in the 
earth reference system, may be computed from 

<|.!; = 6$.sin2 20(l-e^)e^ (5) 

where = 8.56 x 10^*^ cm^ s-^ is the solar neutrino flux at the Earth computed from the Solar 
Standard Model, e is the Monte Carlo branching fraction of visible photons produced in the 
radiative decay (see Fig. Q), 6^ is the neutrino mixing angle, ts^M and Im^e are the average 
flight times of solar neutrinos (assuming an average energy of about 325 eV) from the Sun to 
the Moon and from the Moon to the Earth, respectively. 

The obtained preliminary lifetime lower limits on a possible radiative decay of solar neutrinos 
are given in Table |]. The limits are given as function of mi, Am^ and sin^ 26. 

The data did not allow to obtain lifetime limits for all mass-Am^ combinations. This 
suggests an underestimation of the contribution of the Poisson spot and of other background 
sources; thus the lifetime limits in Tab. |l] are conservative lower limits. 

6. Conclusions 

The analysis of the digitised images measured by a video camera used during the total solar 
eclipse of August 11, 1999 gave us the opportunity to test our data reduction software and also 
lead to some interesting results: 

a) we evidenced the presence in our data of the Moon image in the light reflected from the 
Earth; 
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b) after the Moon image removal and the wavelet decomposition of the remaining signal we 
obtained a correlation function (luminosity versus distance from the center of the Moon) 
consistent with Fresnel diffraction (Poisson spot); 

c) after removal of the Poisson spot contribution we obtain preliminary results on the limits 
of the z/2 lifetime, see Tab. |1|. 

Experimental data of higher quality are needed to obtain more stringents limits. 
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